Fractal Dimension in 3d Spin-Foams 
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In this paper we perform the calculation of the spectral dimension of the space-time in 3d quantum 
gravity using the dynamics of the Ponzano-Regge vertex (PR) and its quantum group generalization 
(Turaev-Viro model (TV) [l(|). We realize this considering a very simple decomposition of the 3d 
space-time and introducing a boundary state which selects a classical geometry on the boundary. 
We obtain that the spectral dimension of the space-time runs from ~ 2 to 3, across a « 1.5 phase, 
when the energy of a probe scalar field decreases from high E < Ep to low energy. For the TV 
model the spectral dimension at hight energy increase with the value of the cosmological constant 
A. At low energy the presence of A does not change the spectral dimension. 



Introduction. In past years many approaches to quan- 
tum gravity studied the fractal properties of the quantum 
space-time. In particular in causal dynamical triangula- 
tion (CDT) [l[ and asymptotically safe quantum gravity 
(ASQG) [31, a fractal analysis of the space-time gives a 
two dimensional effective manifold at high energy. In 
both approaches the spectral dimension is I?s = 2 at 
small scales and Vg = A at large scales. The previous 
ideas have been applied in the context of non commu- 
tativity to a quantum sphere and K-Minkowski [sj and in 
Loop Quantum Gravity Q. The spectral dimension has 
been studied also in the cosmology of a Lifshitz universe 
and in Causal Sets [gj]- Spectral analysis is a use- 
ful tool to understand the effective form of the space at 
small and large scales. We believe that the fractal analy- 
sis could be also a useful tool to predict the behaviour of 
the 2-point and n-point functions at small scales [7] and 
to attack the singularity problems of general relativity in 
a full theory of quantum gravity 9| . 

In this paper we apply to the Ponzano-Regge (PR) 
model and to the Turaev-Viro model (TV) [lal, [IH, 
the analysis introduced in [1]. We consider the ap- 
propriate spinfoam model and we use the very simple de- 
cornposition of the 3d space-time introduced by Speziale 
in [l3[. The other ingredient is the general boun dary 
formalism useful to define the boundary geometry [14| . 
All the space-time is approximated by a single tetrahe- 
dron and the boundary state is peaked on the boundary 
geometry of it. 

The paper is organized as follows. In the first section 
we define the framework and we recall the definition of 
spectral dimension in quantum gravity. The analysis in 
this section is general and not strongly related to the PR 
or TV models. The analysis is correct for any spin-foam 
model. In the second section we calculate explicitly the 
spectral dimension for the PR and the TV theories using 
the general boundary formalism to define the 3d quantum 
gravity path integral. 

The Spectral Dimension. The following definition of a 
fractal dimension is borrowed from the theory of diffu- 
sion processes on fractals |l3| and is easily adapted to 
the quantum gravity context. Let us study the diffu- 
sion of a scalar test (probe) particle on a d-dimensional 



classical Euclidean manifold with a fixed smooth metric 
g^v{x). The corresponding heat-kernel Kg{x,x' ;T) giv- 
ing the probability for the particle to diffuse from x' to 
X during the fictitious diffusion time T (this is just a fic- 
titious time and the scalar field in general is just a tool 
to understand the fractal properties of the space-time) 
satisfies the heat equation 

dTKg{x,x'-T) = ^gKg{x,x'-T) (l) 

where Ag denotes the scalar Laplacian: Ag(/) = 
(^-1/2 ,9^(^1/2 ^M"^ ,9^0) The heat-kernel is a matrix el- 
ement of the operator exp(TAg), 

Kg{x,x'-T)^{x'\eMT^g)\x). (2) 
In the random walk picture its trace per unit volume, 

Pg{T) = V-' J d'^x y/^Kgix, X- T) 



y-^Tr exp(rAg) 



(3) 



has the interpretation of an average return probability. 
(Here V = J d'^x y/g denotes the total volume.) It is well 
known that Pg possesses an asymptotic expansion (for 
T ^ 0) of the form Pg(r) = (4^r)-'^/2 ^^^^ ^„ T". 
For an infinite flat space, for instance, it reads Pg{T) — 
(47rr)-'*/2 foj. all T. Thus from the knowledge of the 
function Pg one can recover the dimensionality of the tar- 
get manifold as the T-independent logarithmic derivative 



d\nPg{T) 



(4) 



dlnT 

This formula can also be used for curved spacetimes and 
spacetimes with finite volume V provided that T is not 
taken too large. 

In quantum gravity it is natural to replace Pg{T) by 
its expectation value on a state \'^). Symbolically, 



P(T) := {Pg{T)) 



DgP{T)e 



^S{g) 



(5) 



Given P{T), the spectral dimension of the quantum 
spacetime is defined in analogy with ([?]): 



^dhiP{T) 
d\nT 



(6) 
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We can formally also to replace the equation ([T]) with the 
correspondent expectation value 



dT{Kg{x,x';T)) = {AgKg{x,x';T)). 



(7) 



The Spectral Dimension in Quantum Gravity. In quan- 
tum gravity we define ([5]) the spectral dimension in the 
general boundary formalism. We introduce a gaussian 
state \ipq) peaked on the boundary geometry q = {q,p) 
defined by the metric and the conjugate momentum. We 
can think the boundary geometry to be the boundary of 
a d-dimensional ball. The state is symbolically given by: 



(8) 



The amplitude (O can be defined for a general spin-foam 
model 



{W\Pg{T)\^^) _ j:s„.,W{Sl){Sl\Pg\S2)MS2) 



(W^l*q> 



(9) 



Where W{s) codifies the spin-foam dynamics ^16^. For 
the purpose of the paper we will consider the PR model 
(TV model); the vertex amplitude is encoded in the {6j}- 
symbol, W{s) oc {6j} {W{s) oc {6j}q for TV and q, the 
quantum deformation of the SU{2) group, is related to 
the cosmological constant A by g = exp(2i-\/A^p)- Since 
we are interested in the scaling of the Laplacian to ana- 
lyze the fractal properties of the space-time, we can ap- 
proximate the metric in the Laplacian with the inverse of 
the SU{2) Casimir operator. We recall that in 3d quan- 
tum gravity the Casimir operator is related to the length 
spectrum of a link e in the simplycial decomposition by 
the relation 



Ll = llC\j,) 



lp[je{je 



1) 



(10) 



where the constant is chosen to be c = 1/4 in line with 
. In 3c? gravity we approximate the 3-ball with a single 
tetrahedron and the boundary sphere by the surface 
of the tetrahedron given by the six triangles. We consider 
fixed four of the six representations (j) and we call the 
other two free representations by je (e = 1,2). Follow- 
ing the ideas and notation above we define the operator 
Pg{T) in the following way, 



PjXT) ■■= V^^Tr e 



(11) 



Where Aq is the Laplacian at a lower infrared scale, 
is fixed (for example) to je = ji and 



a. e 



(12) 



The boundary state in the notation above is 

M.^.(j,) = AA-ie-*^=(^'=-^')'+'^^'(^'=+i/2). (13) 

where A/" is a normalization factor. The dihedral angles 
9 = arccos(— 1/3) define the boundary extrinsic geome- 
try for an equilateral tetrahedron. Now we have all the 




FIG. 1: This is the plot of the modulus of the expectation 
value (fT4|) \t] {W\dji (black dots) compared with the ex- 
ponential exp(— TCq |Aoi/[j(j + !)-!- c]) in the case c = 1/4 
(red dots). The expectation value is calculated for T — 1 and 

r = io, o<ji,j2 =s;2j. 




FIG. 

TCi < 40 



2: Plot of the amplitude d} for 1 < j < 6 and 1 < 



ingredients to calculate the expectation value ([9|) using 
(fT2| and (fT3|l . In particular, since the geometry appears 
only in the operator Oe , we can calculate the expectation 
value of this operator. 



ii,2=o 



2j 6 



= li^"^^- + l){6j}e-^bi^"+i)+=l *j(ji,2), (14) 

where we introduced the following notation for the nor- 
malization, ?7^^ := (IFl^tfj). We also replaced the Lapla- 
cian Aq with — |Ao|. Before to calculate the amplitude 
(fn)) we replace Aq with |Ao| oc this assumption 

will be clear later in the paper. 

The result of the calculation (|14p is given in Fig[T]and 
compared with the exponential exp[— TCg | Ao|/(j(j+l) + 
c)] in the case c = 1/4, Cq = 1 and |Ao| cx 1/f. The 
plots in FiglTJare for T = 1 and T = 10. We can observe a 
perfect agreement for j > 4. This agreement is supported 
by the plots in Fig[2]and Fig[3] where the amplitude 
on the left and the function exp(-X|Ao|/[7(i + l)-)-l/4]) 
{X — TCq) on the right coincide for j > 4. In FiglHwe 
plotted a section of (O for j = 6 and X e [0, 100]. This 
section coincides with the function exp{—X\Ao\/[j{j + 
1) + 1/4]) evaluated on j — 6. 

In the range 1 < j < 12 we have interpolated the 
exact result (ITil) numerically and obtained a different 
exponential form of the amplitude. The points data and 



FIG. 3: Plot of the amplitude ([14)1 for 1 < i < 8 and 1 < X < 
40 on the left and of the function exp( - X | Ao | / [ j ( j + 1 ) + 1 /4] ) 



on the right, where X 
agreement for j > 4 . 



TCo This plots show there is good 



FIG. 5: The plot on the left represents the points ob- 
tained from the evaluation of the amplitude ((T?)) for 1 < 
j < 20 (black points). The red points refer to the function 
exp{-T Co Ao/[ j{j + 1) + c]) for Cq = 1 (it is an irrelevant 
constant) and T = 1; it is evident that for small values of the 
representation j the two function are different. The plot on 
the right represents an interpolation of the black points in the 
picture on the left for 1 < J < 12. 



with the scale £ (or with the energy scale k « 1/^), 



ho A- 



for j>l (j >4), 
for j w 1 il<j< 4). 



(17) 



FIG. 4: Plot of the amplitude ((2} for j = 6 as function 
of X = TCo on the left (black points) and of the function 
exp(-X|Ao|/[i(j + 1) + 1/4]) 1 < TCi < 40. 



the fit are given in FiglHl The points are fitted by the 
function aexp(&/j"), where a « 1.00, b ~ 0.55 and a ~ 
3.03 for T = 1. Recalling that A cx 1/j^ we conclude 
that at the Planck scale, 



(15) 



We will use this result to calculate the spectral dimension 
at the Planck scale then for T « 1 in Planck units, this 
is the reason why we fixed T = 1 in the expectation 
value. We can reproduce the behavior of for j > 
4 (in FiglH the function exp(-TC^Ao/[ j(i + 1) + c]) 
coincides perfectly with the exact expectation value ITU) 
from j > 4) also analytically using the asymptotic large 
j limit of the {6j} symbol. For large j we have: {6j} oc 
exp{iSii{je) + in /A) + c.c. Using this property of the 
symbol and replacing the sum in (|14p with an integral on 
Sji,2 ji,2 ~i (for j > 1) we obtain exp(rC^Ao/[i(j + 
l) + c]). 

What we learnt from the explicit calculation of ([9]) can 
be summarized as follows, 



(a 



e^^iO+iVc^" for j > 1 (j > 4), 
a e-^'-p^"^" for j « 1 (1 < j < 4). 



(16) 



Where a « 1.03. We introduce a DifF-invariant scale 
defined hy £ :— jelp The result of the ^ can be summa- 
rized in the scaling property of the Laplacian operator 



Where we introduced the infrared scales Jo 7 Jo io ^ 
1, Jo ^ 4), c' = bC§ and, by definition, C]^ = jo(jo + 1) + 
c. 

We denote the scaling of the Laplacian operator sug- 
gested by (jl7p by a general function in the momentum 
space. We introduce here a physical input to put the 
momentum k in our analysis. If we want to observe the 
space-time with a microscope of resolution I = Ipj (the 
infrared length is £0 '= Ipjo) we must use a (fictitious) 
probing scalar field of momentum k ^ l/l. The scaling 
property of the Laplacian in terms of k can be obtained 
by replacing: I ~ l/Zc, Iq ^ 1/ko and Ip ^ 1/Ep, where 
/cq is an infrared energy cutoff and Ep is the Planck en- 
ergy. We define the covariant Laplacian at the scale k 
introducing the function Sk, 



(18) 



It is straightforward to derive the scaling function from 
(flT)) and using the arguments above 



§fc(/c, fco 




for < fc < Ep. 



(19) 



We added a factor one in the infrared limit to facilitate 
the spectral dimension calculations. The scaling func- 
tion Sk{k,ko) represents also , using the definition of 
the Laplacian, the scaling of the inverse of the metric, 

{gnk=Mk,ko){9nko- 

We suppose that the diffusion process involves (approx- 
imately) only a small interval of scales near k then the 
corresponding heat kernel contains the A^ for this spe- 
cific and fixed value of the momentum scale k. Denoting 
the eigenvalues of — Afc,, by En and the corresponding 
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eigenfunctions by (j)n{x) — {x\En), we have the following 
eigenvalue equation for the Laplacian 

{E„\Ern)= J d^x' y/go{x') 4>l{x')4>n{x') = 5n.n- (20) 

Using (PO)) and the definition ^ we can calculate explic- 
itly the heat kernel Kk{x,x';T) = (a;'|(C'e)|a;). By using 
([TBI). (El), (BHl) and (EOl) we have 

Kkix,x';T) = {x'\e^^^\x) 

= Y,c^l{x')<l,^{x)e-^^^(^ME^. (2f) 

n 

From the knowledge of the propagation ker- 
nel (|2f p we can time-evolve any initial proba- 
bility distribution p{x; 0) according to p(x] T) — 
J d^x' ga{x') K{x, x'; T) p{x'; 0), where 50 the determi- 
nant of l^g^u) ^ ■ If the initial distribution has an eigen- 

function expansion of the form 0) = X]n^"'/'"(^) 
we obtain for arbitrary x, 

p{x-T)= J dV y/goix')K{x,x';T)p{x';0) = 
= ^C„0„(a;)e-^(^^^-«)^"^ (22) 

n 

where we used the wave function normalization (|20p . If 
the C„'s are significantly different from zero only for a 
single eigenvalue En, we are dealing with a single-scale 
problem and then we can identify = En- However, 
in general the C^'s are different from zero over a wide 
range of eigenvalues. In this case we face a multiscale 
problem where different modes 0„ probe the spacetime 
on different length scales. 

If A(fco) is the Laplacian on the flat space, the eigen- 
functions 4>„ = (j^p are plane waves with momentum p'^, 
and they resolve structures on a length scale £ of order 
Hence, in terms of the eigenvalue En = Ep = p^ 
the resolution is £ « l/^/E^. This suggests that when 
the manifold is probed by a mode with eigenvalue En it 
"sees" the metric {g^i,)^^ for the scale k = \fE^. Ac- 
tually the identification k — \JEn is correct also for a 
curved spacetime because the parameter k just identi- 
fies the scale we are probing. Therefore we can conclude 
that under the spectral sum of we must use the scale 
k- = En which depends explicitly on the resolving power 
of the corresponding mode. In eq. P^ . Sfe(fc, k^) can be 
interpreted as S(i?„). Thus we obtain the traced propa- 
gation kernel, 

PiT) = J2 ^ y " = ^"'Tr (eS(-'^'^o)A.o^). (23) 

n {9)fcQ 

It is convenient to choose fco as a macroscopic scale in 
a regime where there are not strong quantum gravity 
effects. 



We assume for a moment that {gtii^}/^^ is an approxi- 
mately flat metric. In this case the trace in eq. (1^^ is 
easily evaluated in a plane wave basis: 

^(^)-/(S^^"''''''- (24) 

where we used the fiat metric {g^i,)ka ~ S^i^ and A^^ \p) = 
-p^\p). 

The dependence from T in (|24|) determines the frac- 
tal dimensionality of spacetime via In the limits 
T — !• 00 and T — > where we are probing very large and 
small distances, respectively, we obtain the dimension- 
alities corresponding to the largest and smallest length 
scales possible. The limits T ^ 00 and T ^ of P{T) 
are determined by the behaviour of §(p) for p ^ and 
p — > 00, respectively. 

The above assumption that {gfj.u} is flat was not neces- 
sary to obtain the spectral dimension at any fixed scale. 
This follows from the fact that even for a curved metric 
the spectral sum can be represented by an Euler- 
Mac Laurin series which always implies ([^5]) as the lead- 
ing term for T — > 0. 

Now we have all the ingredients to calculate the spec- 
tral dimension using ([M)) inside the definition For 
the PR model the scaling function S{p) is obtained from 
(I19p replacing k with p. The spectral dimension for j ^ 4 
or fc < Ep/A increases from Vs « 1.5 to I?s ~ 3 at low en- 
ergy as it is evident from the plot in FiglB] For 1 < J < 4 
oi' Ep ^ ^ ^ Ep/4 using the proper scaling we find 
Vs « 1.98. We conclude that the fractal dimension de- 
crease from the Planck energy to an intermediate scale 
where take the value « 1.5 (for k « Ep/4) and increase 
again to « 3 at low energy (FigE]) . For the TV model we 
have differences only for j < 4 and the result is plotted 
in FiglB] on the right. That plot gives the spectral di- 
mension as a function of the cosmological constant. The 
spectral dimension is in the range 2.00 ^ ^ 2.059 
for 0.001 < A < 0.009 in Planck units. In other words 
the spectral dimension increases with the increase of the 
cosmological constant at the Planck scale. 

Conclusions and Discussion. In this paper we calcu- 
lated explicitly the spectral dimension (2?s) for the 3d 
quantum spacetime using the Ponzano Regge spinfoam 
model. We considered the simplest decomposition of the 
spacetime and we used the general boundary formalism 
to characterize the scaling properties of the expectation 
value for the traced propagation kernel. Using the techni- 
cal simplifications repeatedly used in the graviton propa- 
gator calculations we have evaluated the nonperturbative 
expectation value of the heat kernel. 

In the PR model and for k < Ep/A we have plotted 
Vs as a function of a fictitious diffusion time T or equiv- 
alently as a function of the length scale. We obtained 
three phases: a short scale phase Ip "^l ^ 4lp of spectral 
dimension Vg ~ 2, an intermediate scale phase £ > 41 p 
of spectral dimension Ds = 1-5 and a large scale phase 
with I?s = 3. 
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FIG. 6: The plot on the left represents the spectral dimension 
for i ^ 4 or fc < Ep/i as function of the fictitious time 
T. The dimension at hight energy is 1.5. We have plotted 
T e [0.005, 10^] and used Ej, = 1000, ko = 0.01. The plot 
on the right represents the spectral dimension as function of 
the cosmological constant A in the Turaev-Viro model at the 
Planck scale (1 < j < 9). 



For the TV model the results are equal for k < Ep / 4 
and to feel the effect of the cosmological constant we 
must goes beyond that energy. The spectral dimension 
depends on A as it is evident from the plot in FiglHl 

We interpret the results in the following way. At high 
energy the spectral dimension is I?s < 3 because the man- 
ifold presents holes typical of an atomic structure. The 
cosmological constant basically decreases the number of 
holes increasing the spectral dimension. 
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